AMC 12 atomic mass compilation data extrapolated for atomic masses of nuclei far from the valley of stability

The experimental mass data from the Atomic Mass Compilation – 2012 (AMC12) has been analyzed for two-neutron separation energies (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{S}}}_{{\bf{2}}{\boldsymbol{n}}}$$\end{document}S2n), two-proton separation energies (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{S}}}_{{\bf{2}}{\boldsymbol{p}}}$$\end{document}S2p), double-beta decay energies (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{Q}}}_{2{{\boldsymbol{\beta }}}^{-}}$$\end{document}Q2β−), and four-beta decay energies (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{Q}}}_{4{{\boldsymbol{\beta }}}^{-}}$$\end{document}Q4β−) and plotted against neutron number and mass number, respectively. A new weighted slope method of extrapolation, tested for known and new mass measurements, has been used to obtain the extrapolated mass values with better precision for more than 1100 nuclei far from the valley of stability, out of which more than 100 are being reported for the first time. A comparison has been made with five of the popular mass models with reference to experimental extrapolated masses from the present work and the Atomic Mass Evaluation 2016 (AME16). The extrapolated experimental atomic mass data will be very useful for both experimentalists and mass-model theoreticians, as well as in simulations of astrophysical r-processes.

characterizing the interacting many-fermionic liquid and the model 4,7,9 of many-body theory, which therefore takes into account the shell and the liquid-drop features non-perturbatively. KUTY 10 is a nuclidic mass formula. This is a revised version of the mass formula constructed by Koura et al. and published in 2000 11 . The WS4 12 is the most recent and improved version of the semi-empirical nuclear mass formula based on the macroscopic method together with the Skyrme energy-density functional, the surface diffuseness correction for unstable nuclei, and the radial basis function (RBF) corrections. These global mass models can be regarded as one of the methods for determining the masses of experimentally not (yet) accessible nuclides.
Local interpolations or extrapolations of the experimental mass values have been tried as an alternative and more plausible way to reach the nuclides far from the valley of stability. The Isobaric Multiplet Mass Equation (IMME) is a powerful interpolation/extrapolation method for lighter nuclei. Mass values derived from IMME had been included in older mass evaluations and the recent high-precision mass measurements have sparked renewed interest in this technique 13,14 .
On the basis of the semi-empirical Bethe-Weizsacker 15,16 mass formula, Way and Wood 17 derived graphical presentations of nuclear decay data. Under the assumption that the atomic masses are smoothly varying functions of neutron number (N) and atomic number (Z), as in the semi-empirical mass formula, the lines were extrapolated to unmeasured nuclei. In 1966, Garvey and Kelson presented a relationship based on an independent particle model of the nucleus 18 : This equation relates to six nuclides. With five known masses, a sixth unknown mass can be predicted 19 . These global mass formulas, however, will in general be too approximate to make predictions to a sufficient degree of accuracy.
Audi-Wapstra successfully employed a local extrapolation method 20 based on the systematics and smoothness of the mass surface and its derivatives. These mass measurements and extrapolations were provided as Atomic Mass Evaluations (AME) and have been regularly updated. An interactive graphical program was used to make estimates for unknown masses from the trends in the mass surfaces. However, this procedure encompasses a "subjective" component in the form of individual judgments.
In an attempt to provide extrapolated masses complementary to the extrapolations in AME, we have used an independent data set of mass values, namely the Atomic Mass Compilation 2012 (AMC12) 21 . There are numerous potential applications for the AMC12. Predicting accurate estimations of unknown, imperfectly known, or ambiguous masses; extrapolation toward the drip lines; and testing theoretical models are some of these key applications. Such an attempt has been made in this work. One of our aims in using an independent data set rather than AME mass compilations and a new extrapolation method is that the presence of two complementary, but independent sets of data would boost the confidence in using these extrapolated mass estimates. In addition to this, following the tradition of providing extrapolated masses of unknown masses for nuclides far from stability by the Atomic Mass Evaluations (AME), having published the Atomic Mass Compilation (AMC12), we wanted to provide extrapolated masses based on the AMC12 mass data as an extension of AMC12 work using a different extrapolation method.
We have carried out a point-to-point local extrapolation to systematically take into account the local variation in the trends that are observed in the derivative sheets. The local trend was quantified in terms of the value of the slopes connecting one isotope to the next. The extrapolated slope for the prediction was obtained by averaging the slopes of three of the immediate neighbors. A weighted average of the three slopes was used, with higher weightage placed on the slope closest to the one to be found. The extrapolated slope value was used to find the unknown value of the derivative. This extrapolation method was independently carried out on four derivative sheets. The error in the extrapolated derivative was calculated by considering the error in the value of the derivatives used for the individual extrapolations. Since our extrapolation method is based on systematic local extrapolation, the method is very effective in following local variations in trends present in the derivative sheets. The limitation of the method would be the one where there are conflicting slopes while averaging the three slopes (i.e., absence of trend). As can be observed from the derivative sheets (shown below), consistent trends are observed while plotting individual values. Hence, our method succeeds in accurately predicting the unknown mass values.

Methods
Plotting nuclear masses as a function of N and Z produces a surface in three dimensions. Due to pairing energy, this surface is separated into four sheets that run parallel in all directions 20 . This regularity can also be observed in the derivative sheets. A derivative is any specified difference in the masses of two adjacent nuclei. Remarkably smooth trends can be observed in the derivate sheets and this regularity can be used for obtaining unknown masses of nuclides by extrapolation from the well-known mass values on the same derivative sheet. The derivative sheets have the advantage of displaying much smaller variations.
Analysis of experimental Atomic Mass data from AMC12. The simple trends in the mass surfaces have been used to obtain unknown masses by observing the regularities in derivatives mentioned below. We consider four different mass sheets: S n 2 as a function of neutron number, S p 2 as a function of proton number, Q 2β − as a function of mass number and, β − Q 4 as a function of mass number. Separation energies in MeV of particles or groups are obtained as the following combination of atomic masses from the AMC12 atomic mass data: Two representative mass surfaces are shown in Figs. 1 and 2, clearly showing the trends in the mass surface. The mass surfaces are calculated using the experimental mass values from the AMC12 mass compilation before extrapolation. Figure 1 shows a S 2n mass sheet in which S 2n values are plotted as a function of neutron number from 60 to 85. Figure 2 shows a Q 2β mass sheet in which Q 2β values are plotted as a function of mass number from 32 to 65. The derivative values prior to extrapolation are represented by the filled circles in both graphs. These graphs were used for extrapolation of unknown masses. Precise estimates of unknown masses of nuclides far from the valley of stability have been obtained by extrapolation from well-known mass values on the same sheet.
For example, the two-neutron separation energies, S n 2 (Z, A) have been computed from the ground state nuclear masses M(Z, A) and M(Z, A-2) from the AMC12 and the neutron mass. The evolution of S n 2 as a function of neutron number shows the well-known regularities (see Fig. 1). For a given proton number, S n 2 decreases smoothly as the neutron number increases. The curves of various isotopic chains are roughly parallel to each other. Similar calculations for S p 2 , β − Q 2 , and β − Q 4 were performed using the ground state atomic mass data from AMC12 21 for the respective nuclides. The S p 2 energies versus neutron number, Q 2β − energies versus mass number and, Q 4β − energies versus mass number have been plotted. The points of isotopic chains have been joined to show the roughly parallel nature of the lines in all cases. All the four sets of plots have been subjected to the extrapolation procedures described here under. In the plots, the experimental data calculated from AMC12 is shown as darkened circles, while the extrapolated data is shown as open circles. The plots in the case of It is quite interesting to see the beauty and the consistency of the roughly parallel behavior even with respect to the extrapolated data points, thus showing the quality of the extrapolation method followed in the present work.
In order to estimate the unknown mass values, extrapolations were performed by independently following the trends in the four derivative sheets mentioned above. The value of the extrapolated derivative was then used to obtain the previously unknown mass value. It was observed that point-to-point extrapolation has to be taken into account in order to meticulously account for the local trends in the derivative sheets. In each of the sheets, trends from the three previous (or following) derivatives were used to estimate the value of the unknown derivative. To further improve the accuracy of the extrapolated derivative, it was ensured that the local trends www.nature.com/scientificdata www.nature.com/scientificdata/ were closely tracked by assigning a higher weightage on the trends from the closest neighbors to the derivative that was to be found. The details of the method used to find the extrapolated derivate and the error in the corresponding extrapolated mass value are provided here. The computational technique developed to obtain the extrapolated derivative is illustrated below with the S 2n derivatives as an example.
Extrapolation: The weighted slope method of computing the extrapolated S 2n values. S n 2 i values in Table 1 represent the two-neutron separation energies, where we want to extrapolate to find the unknown two-neutron separation energy labeled S n 2 10 . For a given nucleus, the table is listed as a function of the neutron number (N i ) and proton number (Z j ). To obtain the extrapolated S 2n value (S n 2 20 in this case) for a nuclide, we use the trend in the sheet by using the existing derivatives from the already known elements. We have considered point-to-point local extrapolation, in order to follow the trends in the derivative sheets very closely. In our method to obtain the extrapolated S 2n value for a nuclide, we used the trend in the sheet by using the already known slopes of the line joining the S 2n values for the next three nuclides.
We use the trend in the S 2n sheet by finding the slopes of the line joining the S 2n values for the neutron numbers N 2 and N 1 of (in this case) the next three elements. The slopes for these three known sets of values are calculated as  23 13 In order to follow the trend with improved accuracy, we have provided more weightage for the slope calculated from the S 2n values immediately next to the one we want to find, than for the slope of the others. As a result, the slope m was found as

Fig. 2 A representative Q 2β
mass-sheet before extrapolation. www.nature.com/scientificdata www.nature.com/scientificdata/ Higher weightage is provided for m 1 so that trends in slopes closest to the nuclei have more influence than if the three slopes were provided with equal weightage. This makes sure that the local variation in the trend lines is followed closely. The c-intercept for the extrapolated value was obtained by substituting the slope, m to the already known S 2n value for the nuclide which is, S n 2 20 and rearranging the following equation: to get the c-intercept: Thus, we obtain the extrapolated S 2n value as, A similar method has been extended to find the unknown S 2n value, depending on where in the sheet of derivatives the method has to be used. In all, this method of local extrapolation by weighted slopes has been used in three different forms, in addition to the way in which it is used above. The second case is similar to the one above, but the two-neutron separation energy that has to be found by extrapolation is S n 2 23 , with all the other values in the above table known. For this, we again use the same weightage for slope, where the weightage progressively decreases for the slope that is away from the proton number for which the two-neutron separation energy has to be found. Once the weighted slope is found, the c-intercept is found using the known two-neutron separation energies of the next neutron number (S n 2 13 in this case) and the slope. From the calculated slope and the c-intercept, we then find S n 2 23 . The other two cases are the ones where, for a given neutron number, the two-neutron separation energies of the nuclei with a proton number one higher and one lower are already known. In this case, we perform interpolation to find the two-neutron separation energy using a calculated slope. We first calculate the slope m 1 by using the two-neutron separation energies with a proton number one lower and with a neutron number the same as that of the nuclei for which the two-neutron separation energy has to be found, and the neutron number one higher. Similarly, we find the slope m 2 from the two nuclei with proton number one higher and the neutron number one lower, and with a neutron number the same as that of the nuclei for which the two-neutron separation energy has to be found. Then the c-intercept is found by using the average of these slopes and the two-neutron separation energy of the nuclei with the same proton number as that for the nuclei for which the two-neutron separation energy has to be found but with a neutron number one lower. Thus, we can interpolate to find the unknown two-neutron separation energy by knowing the slope and the c-intercept as before. The last case is the one that is similar to the one described above, but we also know the two-neutron separation energies of the nuclei with a neutron number lower than the nuclei for which the two-neutron separation energy has to be found and the proton number is lower than, same as, and higher than the nuclide for which the two-neutron separation energy has to be found. That is, we know the two-neutron separation energy of all the nuclei that are around the one whose value we have to interpolate. For this case, we not only find the two-neutron separation energy as we did for the previous case but also find the two-neutron separation energy by using the known values with a neutron number that is the same as, and one lower than, the nuclei for which the two-neutron separation energy has to be found. Then we take the average of the two interpolated two-neutron separation energies.
Similarly, we use these methods to extrapolate (and interpolate) the corresponding values in the S 2p , Q 2β −, and Q 4β − sheets.
Calculation of errors. The error in mass is a very important quantity to know the validity of the mass value that has been determined by extrapolation. The value of the extrapolated derivative relies on the previously known mass value from which the known derivatives were calculated. The error in the previously known mass value is used to obtain the error in the extrapolated derivatives. The steps involved in determining the error in the extrapolated mass value are described below.
The error in the extrapolated S 2n , for example, S n 2 10 is calculated from the method used for extrapolation, where we had The error in m, the slope, Δm is calculated using errors in already existing S 2n values, used for the calculation of the slope:  10 20 Where the ΔS n 2 20 is the error in the known derivative value. Then, the error in the extrapolated mass is given by Where the subscript e in ΔM e (A, Z) denotes that it is the error in the extrapolated mass value. Similar calculations have been carried out with different methods of extrapolation, as described above, for each individual nuclide. The procedure has been then, applied for S 2p , Q 2β − and Q 4β − extrapolated data. The extrapolated data of S 2n , S 2p , β − Q 2 , and Q 4β − were used to calculate the respective extrapolated mass values and thus mass excess values, as well as the uncertainties associated with them. We calculated the weighted average of all the available extrapolated mass excess values and the corresponding error for each nuclide using the standard formula for calculating the weighted averages. The data has been arranged in a tabular form in data set 2 22 , with mass number A, the name of the nuclide, the corresponding atomic number Z, and neutron number N in four consecutive columns. This is followed by the weighted average of extrapolated mass excess value and error designated as AMC12-ME, Error in the next two columns. The next two columns contain the extrapolated mass excess values and error corresponding to that from the AME16 mass evaluation data 23 designated as AME16-ME and Error.
In order to provide a comprehensive comparison, we have collected the theoretical mass model values from four popular mass models: the Duflo-Zuker model (DU-ZU) 6 , the finite-range droplet model (FRDM12) 4,5 , the Hartree-Fock-Bogoliubov mean field model (HFB21) 7 , the infinite nuclear matter model (INM) 8 , the nuclidic mass formula composed of a gross term, an even-odd term and a shell term is as a revised version of the mass formula constructed by Koura et al. (KUTY) 10 and published in 2000 11 , and the most recent WS4 model by Wang et al. 12 . These values corresponding to each nuclide are arranged in the next six columns in Table 2, designated as DU-ZU, HFB21, FRDM12, INM, KUTY, and WS4 respectively.

Data Records
All the data generated in the study is available at Figshare 22 . In the present work, more than 1100 experimental data-based mass predictions have been made, in most cases with reduced uncertainties when compared to AME16 data. In the majority of cases, our current data agrees with AME16 data, but with greater precision. More than 370 new masses near the driplines with relatively low uncertainties where there has been no experimental mass data are being presented for the first time. We have mass predictions for more than 100 nuclides in the 1 ≤ Z ≤ 15 region where there were only 50 reported values from AME16 data. The majority of the new data comes from highly neutron deficient nuclei in the lower mass range and highly neutron rich nuclides in the heavy mass range.
In data set 2 22 , we have presented our extrapolation data along with the only available exhaustive extrapolation data from AME16. As there is no way of testing our predictions, we tried to test for self-consistency in data set 3 22 , and in data set 4 22 predictability by using the most recent experimental mass measurements. In data set 5 22 , we have considered a typical region of As, Se, Br, Kr, Rb, Sr, and Y to compare our results with two of the available extrapolations. In the same table, we have also compared our data with one of the most recent mass models, WS4. Four of the five considered models have been widely used for a relatively long time, but the INM model is relatively new and has not yet been tested for its predictive power. Thus, the present work is a real test of the predictive power of these models.

DU-ZU
Mass excess values from the Duflo-Zuker mass model whereby assuming the existence of a pseudo potential smooth enough to do Hartree-Fock variations and good enough to describe nuclear structure, they construct mass formulas that rely on general scaling arguments and on a schematic reading of shell model calculations

FRDM12
Mass excess values from the Finite Range Droplet Model where the calculations are based on the finite-range droplet macroscopic and the folded-Yukawa single-particle microscopic nuclear-structure models, which are completely specified the with considerably improved treatment of deformation and fewer approximations.

INM
Mass excess values from the infinite nuclear matter model, in which the classical liquid drop is replaced by an INM sphere defining the interacting many-fermionic liquid.

KUTY
Mass excess values from the nuclidic mass formula composed of a gross term, an even-odd term and a shell term is as a revised version of the mass formula constructed by Koura et al.

WS4
Mass excess values from the WS4 model, by Wang et al., an improved version of the macroscopic-microscopic mass formula together with the surface diffuseness correction for unstable nuclei, plus radial basis function (RBF) corrections are combined in the WS4 calculations.
www.nature.com/scientificdata www.nature.com/scientificdata/ Data set 2, titled "comparison of extrapolated mass excess values from AMC12, AME16, and mass model predictions", is available through Figshare 22 . The starred (*) M.E. values in the AME16 extrapolated column are the newly measured data reported by AME16 evaluators and are also used in our present work for confirmation of the efficacy of the present extrapolation method in data set 4.

Explanation for data set 3
Data set 3, titled "comparison of extrapolated S 2n values from AMC12 with the known S 2n values" is available through

Technical Validation
In order to test the efficacy of the present method, we tried to extrapolate for the S 2n values of a few known nuclides and compared them with the actual derived values from AMC12 and shown in Table 3. It can be seen that our extrapolation method accurately reproduces the original data. Also, we tried to find how our present extrapolations work in terms of predictability by comparing them with the most recent 24

∆ME
The difference between mass excess data.

Present
The difference between the present extrapolation and the new measurement.

AME12
The difference between the AME16 extrapolation and the new measurement.

WS4
The difference between the value from WS4 and the new measurement. ∆ PW-WS4 Deviation of the present work from the WS4 mass model predictions.
In order to thoroughly evaluate our model, we have run an example test by specifically taking only a subset of known values for As, Se, Br, Kr, Rb, Sr, and Y nuclides for which predictions exist from Garvey and Kelson's global mass formula and the AME16 extrapolations, as well as the most recent WS4 mass model predictions. In Data set 5, we compare our model predictions with these three predictions. Our predictions very well agree with all three extrapolations within our uncertainties.
While the Isobaric Multiplet Mass Equation (IMME) is a powerful interpolation/extrapolation instrument only in the case of lighter nuclei, Garvey and Kelson's global mass formulas will in general be too approximate to make predictions to a sufficient degree of accuracy, and the interactive graphical procedure developed by Audi and Wapstra encompasses a "subjective" component in the form of individual judgments. Even though we don't claim superiority of our present local extrapolation method, it is emphasized that the present local extrapolation method does not suffer from the limitations of the earlier methods. From the success of our extrapolation method in predicting the unknown mass values, which can be seen from the comparison with other models, the present method of using weighted slopes to compute S 2n values can be said to be an effective method to obtain dependable estimates of unknown, poorly known, or questionable masses, and it also does not depend on any personal judgment.

Usage Notes
An attempt has been made to improve the efficacy of local extrapolations of experimental mass data through a simple weighted slope method. The predictive power is remarkably encouraging, as can be seen from the comparison with very recent new mass measurements, where the deviations are well within the experimental uncertainties in most cases. Our extrapolations, in the majority of cases, show very good improvement over the earlier extrapolations.
As can be seen from data set 2 22 , where we have presented our experimental extrapolation data with the only available exhaustive extrapolation data of AME16, our method could predict further away from stability in some cases, and for the majority of the nuclei, the accuracy of our data is either better or comparable. The self-consistency of our method of extrapolation is well tested in data set 3 22 and in data set 4 22 , where we compared our data with the most recent experimental results available, demonstrating the predictability well within the experimental uncertainties. A humble attempt is made in data set 5 22 , where we have considered a typical region of As, Se, Br, Kr, Rb, Sr, and Y and compared our predictions with two available extrapolations and one of the most recent mass models, WS4, we can see the deviations are well within the uncertainties. In our data set 2 22 we have included our extrapolated data, AME16 extrapolated data and the mass model predictions for each nuclide so that the experimentalists can have the full information at hand in one place. This would provide an opportunity for the theoreticians to look for new systematics and trends in different regions of interest and accordingly improve the calculations, while the experimentalists can plan new measurements. Interestingly, with our limited analysis of the comparison with mass model predictions, the most recent WS4 calculations together with the surface diffuseness correction for unstable nuclei and the RBF corrections show very good agreement within the experimental measurements, better than the HFB and Duflo-Zuker formulas. Thus, the extrapolated experimental data of masses can be useful for assessing the impact of current and future experiments in the context of model developments. Additionally, it is anticipated to have an impact on research in fields where experiments are currently not viable, including in simulations of the astrophysical r-process.

Code availability
The extrapolation can be performed by a straightforward implementation of the equations provided in the methods section. To aid the implementation of the equations, an example extrapolation is provided in spreadsheet (.xlsx) format in Figshare 22 . The illustration depicts the extrapolation method to obtain the S n 2 10 value in Table 1, and the extrapolated S 2n value can then be used to obtain the corresponding unknown mass value.